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ABSTRACT 

We show that, unhke the results presented previously in the literature, 
the transition from an outer Shakura-Sunyaev disk (SSD) to an advection- 
dominated accretion flow (ADAF) is possible for large values of the viscosity 
parameter a > 0.5. The transition is triggered by thermal instability of a 
radiation- pressure-supported SSD. The transition radius is close to the central 
black hole. We confirm our qualitative prediction by actually constructing 
global bimodal SSD- ADAF solutions. 

Subject headings: accretion, accretion disks - black hole physics - hydrodynamics 
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1. Introduction 

The model of bimodal accretion disk, which consists of a geometrically thin, optically 
thick Shakura-Sunyaev disk (SSD) (Shakura & Sunyaev 1973) as the outer part and a 
quasi-spherical, optically thin ad vection- dominated accretion flow (ADAF) as the inner 
part, has been quite successfully applied to black hole X-ray binaries and galactic nuclei 
(see Narayan, Mahadevan, & Quataert 1998 for a review). In this model, the accretion flow 
switches from an SSD to an ADAF at a transition radius Rtr- However, Rfr was infered only 
by making plausible assumptions on how it depends on the accretion rate M. The precise 
mechanism through which the SSD material is converted into an ADAF remains a matter 
of debate, and no self-consistent global solution for such a transition has been found so far. 
Indeed, Narayan, Kato, & Honma (1997) obtained a number of examples of global transonic 
ADAF solutions which connect outward to geometrically thin disks. Lu, Gu, & Yuan (1999) 
recovered the whole class of such global ADAF-thin disk solutions. But these authors did 
not take the local radiative cooling into account, thus their solutions could not show the 
variation of the optical depth and could not be regarded as bimodal SSD- ADAF solutions. 
On the other hand, Chen, Abramowicz, & Lasota (1997) considered the local radiative 
coohng as provided by thermal bremsstrahlung, but did not find SSD- ADAF solutions. 
Igumenshchev, Abramowicz, & Novikov (1998) also adopted bremsstrahlung cooling and 
found that only an outer Shapiro-Lightman-Eardley (SLE) disk (Shapiro, Lightman & 
Eardley 1976) could smoothly match an inner ADAF. Such SLE-ADAF solutions are 
optically thin everywhere. More impressive in this respect is the result of DuUemond & 
TuroUa (1998, hereafter DT98). They concluded that the SSD- ADAF transition was not 
permitted, and only an outer SLE disk could match an inner ADAF or SSD. But, as they 
noticed, because of the thermal instabihty of the SLE disk, such SLE-ADAF or SLE-SSD 
models most probably do not exist in nature. 
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It would be a pity if such a promising bimodal SSD-ADAF model could not be actually 
constructed. In this Letter we also use an argument based on energetic considerations 
which is similar to DT98, but find that the SSD-ADAF transition is possible for large values 
of viscosity parameter. We confirm our qualitative prediction by presenting a numerical 
example of the global SSD-ADAF bimodal solution. 



2. Equations 

The system of equations adopted here, which is similar to DT98, makes the following 
standard assumptions: 

(1) : The vertical half-thickness of the disk is expressed as H — Cs/O^k, where Cg — {p/ pY^"^ 
is the isothermal sound speed, with p and p being the total pressure and the mass density 
at the equatorial plane, respectively, and Vt^ is the Keplerian angular velocity calculated by 
using the Newtonian potential, VLk — {GM/R^y/'^. 

(2) : The kinematic viscosity coefficient is expressed as u — aCgH. 

(3) : p is the sum of gas and radiation pressure, p = Pg + Pr- Pg = pT^iTi + Tg), where 
Tj and Tg are the ion temperature and the electron temperature, respectively, and 

Tg =min(T'j, 6 x l^fK). pj. — Qradi'T + 2/-\/3)/4c, where Qrad is the radiative cooling rate 
and T = KpH is the total optical depth. 

(4) : The opacity n is the sum of electron scattering and absorption opacity, k = Kes + Habs, 
where Kgs = OMcm^g~^ and Kabs = 0.27 x 10"^^ pT'^-^cm^g'^. 

Introducing the radial velocity vr, the angular velocity Q, and the surface density 
E = 2i?p, the continuity, radial momentum, azimuthal momentum, and energy equations 
take the form (e.g., Narayan & Yi 1994; DT98) 

M = -2nJ:RvR , (1) 
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i*.0. (2) 



— Q — Qadv + Qrad ■ (4) 



The viscous heating Q'^ has the usual expression 



The advective coohng is expressed as 

n \- r^ds 1 del dp, 

Qad. = ^vnT- = ^vn{^^^ - --) . (6) 

The radiative coohng rate Qrad is calculated using a bridging formula (DT98), 

aa. = 8aT,^(^ + V3 + ^)-. (7) 

Eq.(7) is vahd in both optically thin and optically thick regimes. The bremsstrahlung 
coohng is given by (Abramowicz et al. 1995) 

Q^em = 1-24 X IQ^^Hp^Tl/'^ergs s-^cm~'^ . (8) 



Abramowicz et al. (1995) and Chen et al. (1995, hereafter Chen95) obtained a unified 
S — M picture of accretion flows around black holes with an additional assumption that 
the disk always rotates at O = Q.k- However, it is known that slim disks (Abramowicz 
et al. 1988) and ADAFs are quite sub-Keplerianly rotating. We therefore assume: 
Q = ujVLk{^ < < 1), where the parameter uj is obtained from the self-similar solution 
(Narayan & Yi, 1994). We fix 7 = 1.5 and adopt the self-similar method, then the equations 
(2), (3), (5) and (6) are reduced to the following algebraic form: 

\vl + + {u^ - mlR^ = , (9) 



Q+ = Am02 ^ 



(11) 



o - 1 no) 

Q^v - . (12) 

Prom Eq.(9) one easily obtains: h = H/R = Cs/flxR < ^2/5 fa 0.63. We use a 
parameter / = Qadv/Q^ = /'2>uj'^ to measure the degree to which the flow is advective. If 
the radial velocity is neglected, Eq.(9) is simplified to: (o;^ — 1) + — 0. For radiative 
coohng-dominated flows (/ ^ 0) such as SSDs or SLE disks, u ^ 1, and /i 0; while 



for advection-dominated flows (/ ~ 1) such as slim disks or ADAFs, oj ~ w2/17, and 



h ~ ^J^Jvi . Thus the assumption ^ = ujVLk is more reasonable as it is valid for both 
radiative cooling-dominated and advection-dominated accretion flows. 



3. SSD-ADAF Solutions 

Abramowicz et al. (1995) flrst obtained a unifled E — M picture for accretion flows at 
a flxed radius R in the case of low viscosity, which includes four classes of solutions, namely 
SSDs, SLE disks, slim disks and ADAFs. Chen95 found that two types of E — M picture 
should exist, which are separated by a critical viscosity parameter acrit- We recover here 
these two types of E — M picture with our assumption Q — cuD,k, which are given in Fig.l, 
(a) for a < acrit and (b) for a > acrit- The solid lines in the figures represent thermal 
equilibrium solutions, i.e. with — Q~ . In Fig. 1(a), the right ^'-shaped curve consists 
of three branches, of which the lower one is for gas-pressure-supported SSDs, the middle 
one for radiation-prcssurc-supported SSDs, and the upper one for slim disks; while the left 
curve consists of two branches, of which the lower one is for SLE disks, and the upper one 
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for ADAFs. In Fig. 1(b), the straight hne is for ADAFs and shm disks, while the n-shaped 
curve consists of three branches, of which the two branches on the right are the same as the 
middle and lower branches of the ^-shaped curve in Fig. 1(a), and the branch on the left is 
for SLE disks. The unstable branches are those which have > Q' above and Q'^ < Q~ 
below, while the stable branches are just opposite. Thus gas-pressure-supported SSDs, slim 
disks and ADAFs are thermally stable, whereas radiation-pressure-supported SSDs and 
SLE disks are thermally unstable. 

DT98 investigated the thermal instability in Fig. 1(a), and argued that the SSD-ADAF 
transition is not permitted. We agree with DT98 that such a transition is indeed prohibited 
in Fig. 1(a). The arrows in Fig. 1(a) show a hmit-cycle behavior resulting from the thermal 
instability of a radiation-pressure-supported SSD. 

However, the other type of picture as shown in Fig.l(b) was ignored by DT98. Because 
a bridging formula like Eq.(7) is used, the optically thick, high M slim disk solution and 
the optically thin, low M ADAF solution can be described by a single line in Fig. 1(b), the 
line extends over the entire range of M without break. It is this feature that makes the 
SSD-ADAF transition possible. A thermal disturbance on a radiation-pressure-supported 
SSD can trigger the flow to behavior in the following way: The flow flrst jumps to a slim 
disk solution and becomes thermally stable. But because the accretion rate M does not 
match that of the outer SSD, then the shm disk must evolve into an ADAF, for which 
M matches that of the outer SSD. The whole process is indicated by the two arrows in 
Fig. 1(b). An SSD-ADAF transition is realized. 

To confirm this qualitative prediction, we now go on to present our numerical models 
for bimodal SSD-ADAF disks. The general thermal instability condition is 



We denote (3 — Pg/ipg +Pr), and A = K,abs/{>^es + >^abs)- In privious researches on SSDs the 




(13) 
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disk was usually divided into three separate regions, for which (from the inner to the outer) 
the two parameters A) are (0,0), (1,0) and (1,1), respectively. Here we let both j3 and 
A vary continuously from to 1, thus the SSD solution obtained will smoothly extend over 
these three regions. In our formulation Eq.(13) takes an explicit form 

(5 = 4 - 10/? - 7.5A - 0.5/3A > . (14) 

When the opacity is dominated by the electron scattering, i.e. A = 0, this condition is 
reduced to the usual expression: (3 < 0.4. 

Once the viscosity parameter a and the relative accretion rate m {fn — M/MEdd 
with MEdd being the Eddington accretion rate) are given, we can search for the global 
SSD-ADAF solution by the following steps: 

(i) Provided that the given viscosity parameter a > 0.5, a critical radius Rc exists. The 
E — M picture is of the type of Fig.l(a) for R > R^ and of the type of Fig. 1(b) for it! < Rc- 
We first calculate the value of Rc which corresponds to the given a. 

(ii) We solve the equations of SSD inward from an outer boundary Rout = 2000Rg. The 
SSD solution breaks off when the thermal instability condition Eq.(14) is met at a certain 
radius Rf,. We calculate the value of Ri, which corresponds to the given a and m. If the 
condition Eq.(14) is never met for the given a and m, then the SSD solution is thermally 
stable everywhere. 

(iii) If Rb < Rc, the SSD-ADAF transition will occur at Rtr — Rb- We can obtain the 
self-similar ADAF solution inside Rtr. Combining with the SSD solution outside Rtr, we 
obtain the whole global SSD-ADAF solution. On the contrary, if Rh > Rc, the limit-cycle 
behavior will occur at R^, and the SSD-ADAF transition will not occur. 

In Fig. 2 we show how the behavior of an original SSD flow depends on a and fn. It 
is seen that three possible cases, namely the stable SSD, the limit cycle behavior, and 
the SSD-ADAF transition, each occupy a certain region of the a — m parameter space. 



-9- 



and the SSD-ADAF transition occurs for a > 0.5 and a definite range of m. To construct 
numerically a global bimodal SSD-ADAF solution, we choose a = 0.8. We then obtain 
the corresponding critical radius Rc = 12.0Rg, and find from Fig. 2 that the SSD-ADAF 
transition occurs for 0.03 < m < 0.17. The SSD is stable everywhere for m < 0.03, and the 
limit-cycle behavior appears for m > 0.17. An example of the global SSD-ADAF solution 
for m = 0.1 is presented in Fig. 3. The SSD becomes unstable at Rb = S.ORg. The figure 
gives a logijl) — log{rh) picture aX, R — 8.0Rg, which belongs to the type of Fig.l(b). The 
arrows indicate the transition from the SSD solution (filled circle) to an ADAF solution 
(filled square). 



4. Discussion 



In this Letter we show that the thermal instability of a radiation-pressure-supported 
SSD can possibly trigger two different kinds of behavior of the flow, namely the limit-cycle 
and the SSD-ADAF transition. For low values of viscosity parameter a < 0.5, only the 
limit-cycle behavior can occur; while for large values a > 0.5, either of the two kinds 
of behavior can, and which one is actually realized is determined by rh. We use two 
parameter, /3 and A, to smoothly connect the three usually separated regions of an SSD, so 
that the exact position of the SSD-ADAF transition can be found, and the global bimodal 
SSD-ADAF solution can be obtained. 

The range of a values necessary for the SSD-ADAF transition to occur should depend 
on, among other things, radiative cooling mechanisms adopted. Chen95 used two bridging 
formulae of Qrad to calculate acru for a flxed radius R, one of which is introduced by 
Narayan & Yi (1995), who considered bremsstrahlung and synchrotron cooling, and 
Comptonization; and the other of which is introduced by Wandel & Liang (1991), who 
considered bremsstrahlung cooling with Comptonization only. They found that the value 
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of acrit calculated with the former formula (i.e. with more sources of radiative cooling) 
is higher than that calculated with the latter one. For example, for R = 5Rg, they had 
Q^crit = 0.41 with the former formula, and acrit = 0.245 with the latter one. Since in bridging 
formula Eq.(7) only bremsstrahlung cooling is taken into account, we guess that if other 
sources of radiative cooling such as synchrotron emission and Comptonization are included, 
the lowest required value of a for the SSD-ADAF transition would become larger than 0.5, 
i.e. the necessary condition on a for the existence of bimodal SSD-ADAF disk would seem 
even more restrictive. It is unclear yet if there are other factors which could help to relax 
this condition on a. However, the good thing for the SSD-ADAF transition is that the 
larger the allowed value of a is, the wider the corresponding range of rh is, as seen clearly 
from Fig. 2. 

The SSD-ADAF transition found in this Letter results from the instability of the 
radiation-pressure-supported region, i.e. the inner part of SSD. Thus the transition 
position is close to the central black hole, as we have calculated here. On the other 
hand, observations of some black hole X-ray binaries (Narayan et al. 1998 and references 
therein) seem to imply a SSD-ADAF transition radius Rtr ~ iO'^Rg, which is definitely in 
the very outer region of SSD. However, it is well known that the outer region of SSD is 
gas-pressure-supported and is both thermally and viscously stable. In our opinion, it is 
difficult to see from the theoretical point of view how to convert a gas-pressure-supported 
SSD directly into an ADAF. The apparent conflict between present observations and 
theories remains an unsolved issue. 
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Fig. 1. — Two types of logCE) — log{M) picture for accretion flows at a fixed radius. The 
Solid lines mark the different branches of solutions. The arrows show the behavior of the flow 
resulting from the thermal instability of a radiation-pressure-supported SSD. (a) a < acrit, 
a hmit-cycle occurs, (b) a > acrit, an SSD-ADAF transition occurs. 

Fig. 2. — Dependence of the possible behavior of an original SSD flow on a and m. 

Fig. 3. — An example of global SSD-ADAF solution with a = 0.8, and m = 0.1. The 
logiTi) — log{m) picture is for the accretion flow R — S.ORg, where the SSD becomes 
thermally unstable. The arrows show the transition from the SSD (filled circle) to an ADAF 
(filled square). The dotted line corresponds to the limit: H/R = 0.63. 



